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• Abstract. We study two-body correlations in a many-boson system with a 

j. . hyperspiierical approach, where we can use arbitrary scattering length and 

O ! include two-body bound states. As a special application we look on Bose- 



o 



Einstein condensation and calculate the stability criterium in a comparison 
with the experimental criterium and the theoretical criterium from the Gross- 
ed ! Pitaevskii equation. 

o 

^ 1 Introduction 

^ ■ Bose-Einstein condensates (BEC) of alkali atoms, e.g., Rb and Na, have suc- 

O ■ cesfully been created in the last decade The JILA experiments jH E] touch 

upon the collapse of boson systems at large two-body scattering length due 
to the formation of two-body bound states. The stability criterium has recently 
^ . been measured to be A^|ac,|/6o < 0.55 |4j, where is the number of bosons, Ug is 

^ I the two-body s-wave scattering length, and 6q = y^h/(mu}) with the geometric 

mean cu of the angular frequencies of the deformed oscillator. The properties of 
BEC have largely been accounted for by the Gross-Pitaevskii equation (GPE) 
[31 E] with a zero-range interaction, which reproduces the correct average energy, 
density distribution, and stability criterium This approximation is usually 
succesful at low density n, where the combination nla^p is very small. 

The need for a more realistic potential is becoming more important with the 
cases of stronger attraction or larger densities, where critical phenomena like 
collapse of a condensate j2j and the conversion into molecular BEC occur [2]. 
In addition, the pure mean-field structure of the wave function is not able to 
account for the correlations between the particles. Attempts on top of the mean- 
field have been tried, e.g., including pairing correlations between the bosons 18] in 
the Hartree-Fock-Bogoliubov formulation. Alternatively the Jastrow ansatz with 
the many-body wave function as a product of two-body amplitudes have lead to 
good results for large scattering length j^]. Quantum Monte-Carlo calculations 
have been applied for a large number of bosons in a density matrix formulation 
[TU] . and more detailed studies for smaller numbers PTI IT^ IT^ . 
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Correlated many-boson systems with arbitrary scattering length 



Recently jl4j a hyperspherical description analogous to the mean-field yielded 
the stability criterium A^|as|/6o < 0.67 for a spherical external field. This was 
extended 031 to treat two-body correlations with realistic finite-range potentials 
and applied to a small number of particles. Recently ^H] this was applied to a 
large number of particles and some universal scalings were extracted from the 
results. This paper will review these results and discuss the various structure 
descriptions possible within the model. 



2 Theory 

We study the A^-boson system of identical, interacting bosons of mass m trapped 
by an isotropic harmonic external field of angular frequency uj. The total Hamil- 
tonian is then 

4=1 i=l i<j=l 

We focus on the relative degrees of freedom and use hyperspherical coordinates 
[T^ : hyper radius p given by p'^ = "^i^j rfj / N , where rij = \ri — rj\, and hy- 
perangles O relating to the remaining 3N — 4 relative degrees of freedom. Due 
to the properties of the harmonic oscillator, the center-of-mass motion separates 
out, leaving a relative eigenvalue equation: 

(iJtotal - i^c.m. - E)^{p, Q)=0. (2) 

We factorize the relative wave function as ^{p, J?) = p-^^^'^^^ f{p)^{p, This 
way we obtain an effective radial equation for /(p) as the eigenfunction of energy 
E in the effective potential U (p) 

-^^ + Uip)-E)fip) = 0, (3) 

2m^(p) ^ (3jV-4)(3iV-6) , , ^ ... 

h"^ 4p2 + p2 + ^4 • M 

Here bt = \/h/ (muj) is the length unit of a harmonic trapping potential of 
frequency v = uj/{2tt). The p-dependent A is an effective angular potential, 
which includes the effects of interactions and correlations. 

Assuming then a wave function which is basically a sum, and not a product, 
of two-body amplitudes as 

^ip,n) = Y^<Pip,u,) (5) 

we are able to solve the angular part of the many-body Schrodinger equa- 
tion as just a one-dimensional integro-differential equation with up to two- 
dimensional integrals. In the angular equation we use a finite-range Gaussian 
V{rij) = Vo exp(— rj? /?)^) as the two-body interaction potential. The only ad- 
ditional approximation consists in assuming that the range b of this potential 
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is much smaller than the average distance between the particles. The scatter- 
ing length as can assume any value, and there can be any number of two-body 
bound states, at least in principle; in actual calculations the accuracy is better 
with as few two-body bound states as possible. The usual zero-range interaction 
used in the mean-field amounts to equating with the Born-approximation cb 
to the scattering length. In comparison we generally have that and ob differ 
substantially. The Born-approximation is given by 



When we assume that ^{p, is independent of angular coordinates, the 
angular potential becomes for » 1 



where we due to the mean-field like assumption of no internal structure replaced 
a-Q by Qs- This is generally the asymptotic angular eigenvalue for the lowest 
state above two-body bound states. The mean-field description only contains 
this effective interaction potential. Two-body structures and bound states in 
this two-body correlated model is clearly beyond the mean-field. 

3 Results 

With a Gaussian two-body interaction potential of range h the qualitative fea- 
tures of A are as shown in figure The large-p features only depend on and 
therefore h only matters for the model-dependent region at small p {p/h < 10^ in 
figure^. At small and negative scattering length the angular eigenvalue quickly 
approaches zero as l.bN'^^'^as/ p. At larger, but still negative, scattering length 
the angular eigenvalue follows a constant value — 1.6A^''/'^ before approaching 
zero as before. In the case of a bound two-body state, when the scattering length 
turns positive, A diverges as — Thus a description of two-body bound 
states is also possible within the model. 

The inclusion of the two-body bound state is evident from considering the 
full spectrum of angular eigenvalues for the case of, e.g., one two-body bound 
state. In figure I^J we show the two lowest angular eigenvalues in such a case with 
one two-body bound state and positive scattering length. The lowest eigenvalue 
(dashed curve) diverges to minus infinity proportional to p'^. This corresponds 
to the bound state. The second eigenvalue (solid curve) is negative at small 
hyperradii but turns positive at larger and approaches the asymptotic behavior 
A oc Ug/p (dotted curve, see inset). 

The deviations of this method from the mean- field are illustrated in figure IHl 
where the lowest angular eigenvalue for a case with no two-body bound state and 
negative scattering length is compared to the zero-range angular eigenvalue for 
the same scattering length. For low density n|asp < the effective energy 

of the two methods coincide, for larger densities the mean-field energy diverges, 
while the energy from the finite-range model remains finite. Moreover, it deviates 
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Figure 1. The angular eigenvalue A for A'^ = 1000 and a Gaussian interaction potential of range 
b. The hyperradius p is in units of b and the two-body s-wave scattering lengths as is indicated 
as as/b on the plot. 




in a region where the density is still relatively low n|osp < 1 so higher-order 
correlations (especially three-body) do not play a role yet. 

A direct comparison of the interaction energy is done by calculating the 
interaction energy per particle for the zero-range mean-field and for this finite- 
range model with correlations, see figure |3 The solid curve shows the GPE 
interaction energy, which increases in magnitude until there is no stable system 
for 1 0^1/60 > 0.55. The crosses are the results from this finite-range corre- 
lated model, which almost coincide with the GPE results. Setting the barrier 
height of the radial potential with A = A5 equal to the oscillator energy yields 
the stability criterium A^|as|/6o < 0.53, which agrees well with the above. The 
dashed curve is the result for the mean-field with a scattering length equal to 
the Born-approximation of the finite-range Gaussian we used in the calculation. 
The deviation of the two GPE calculations, is equivalent to the observation that 
the mean- field model with a zero-range interaction gives the correct result, but 
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Figure 3. The lowest eigenvalue for aa/h = —10*, N = 100, and no bound two-body states 
(solid). The dashed curve is for the same scattering length. The vertical lines indicate regions 
of different density. 

not with a finite-range interaction. However, the correlated model reproduces the 
correct interaction energy with a finite-range interaction of the true scattering 
length. We interpret this as a confirmation that the crucial degrees of freedom 
are included in the ansatz for the wave function. 



N\as\/h 




Figure 4. Interaction energy per particle for b^/h = 1442, aa/h = —0.84 (ag/& = —0.5) 
calculated by the GPE (solid line) and the two-body correlated method (crosses). The dashed 
line is the GPE calculation with 0^ = 03 = —0.56. 



4 Conclusion 

We presented the key results of a hyperspherical study of two-body correlations 
in BEC. Wc confirmed the stability critcrium. Moreover, structure beyond the 
mean-field is observed due to the description of two-body bound states and the 
effective energy at larger densities. We expect this model in the future allows 
more detailed analysis of the coupling between molecular BEC and atomic BEC. 



6 Correlated many-boson systems with arbitrary scattering length 

A possible extension is the inclusion of the important three-body correlations for 
a study of the recombination process in dilute atomic medium. 
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